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Abstract 

O 

CN | A random matrix model of black holes is given based on analysis of Gaussian complex ensembles, 
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based on the generalization of chRMT of QCD. Spacetime freedoms are incorporated in terms of 
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eigenvalues of the ensemble. Classical observables of black holes can be computed using this model 
and it satisfies asymptotic freedom and holographic principle. 
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I. INTRODUCTION 



Random matrices are matrices which entries are random variables. Random matrices were 

n n 

first introduced by Hsu [7] and then intensively studied by physicists in nuclear physics [10]. 
Random matrix theory(RMT), has been a very successful theory in the analysis of nuclear 
spectrum. Then it had been applied to condensed matter physics and recently in QCD 
physics, 2D string theory and etc 9] 18]. In {2], the authors conjectured that M theory is 
equivalent to N — > 00 supersymmetric matrix quantum mechanics describing Dq branes. 
RMT was concerned with black hole physics in [13j to estimate the information lost in black 
hole radiation by computing S matrix. In the present work, we shall propose a new idea to 
compute black hole radiation based on RMT. 



II. THE STATISTICAL NATURE OF MEASUREMENT OF TIME 

It has been recognized that there is a deep relation among general relativity, thermo- 
dynamics and quantum mechanics, since classical black hole thermodynamics was found at 
1970s [5J. Although, it is manifest that a full quantum treatment of gravity has to involve 
thermodynamical arguments, the spacetime geometry and thermodynamics are dealt sepa- 
rately in unified theories so far. Since Hawking effect, which is mathematically equivalent to 
Unruh effect, reveals the thermodynamical properties of black holes, which is a new effect 
at large energy scale, compared to classical general relativity. Therefore, thermodynamical 
nature of spacetime or time, has to be considered in a consistent quantum gravity theory. 
ll| is concerned with the fundamental difficulty of quantum gravity from a philosophical 
and historical point of view. 

In the present paper, we would like to relate the nature of time to how we measure time. 
Since 1949, the most accurate measurement of time is based on the mechanism of radiation 
of cesium atoms. One second is defined as the duration of 9,192,631,770 cycles of radiation 
corresponding to the transition between two energy levels of the ground state of the cesium- 
133 atom. Although, this accuracy is far away from plank scale, we still conjecture any 
measurement of time has to involve statistical behavior of a physical system. Therefore 
instead of writing down explicit time coordinate, we consider the statistical behavior of a 
gravitational system when we consider quantum gravity issues. In one word, we believe 
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gravitation is a probabilistic potential. We shall see this physical idea is realized by random 
matrix theory. 



III. SPACETIME FREEDOM IN THE RMT MODEL OF BLACK HOLES 

As we know, at high energies, QCD can be well understood by perturbation theory, by 
the virtue of asymptotic freedom. The chiral symmetry is restored at about 150 MeV by 
lattice QCD stimulation [8]. In the authors showed that constraints imposed by chiral 
symmetry and its spontaneous symmetry breaking determine the structure of low-energy 
effective partition for the Dirac spectrum. The exact results for the low- lying eigenvalues of 
the QCD Dirac operator. The statistical properties of these eigenvalues are universal, which 
is a the property of chiral random matrix model. The partition function of chRMT |l9] is 



N f 

Z^ u ( mi ,... ,m Nf ) = J DWl[det(D + m f )e-^ Tr ^ w+w \ (1) 



where (3 is the Dyson index 



/=i 




>■ " , (2) 



and W + is an n x m matrix with v = n — m and N = n + m. Nf is the number of flavors 
of quarks. 

The masses of quarks can be thought of taking complex values to involve the 9 angle which 
causes CP violation(see 18[]). Therefore D + rrif is neither hermitian, nor anti hermitian 
either. It is tempt to generalize this partition function to black hole states. At extremely 
high energy or high temperature, it is meaningful to enlarge the chiral symmetry group 
SU(Nf) x SU(Nf) to U(2Nf). Use the same W to denote arbitrary 2NfNf complex matrices. 
The partition function of this black hole is 

Z B (N) = J DWdetWe-* Tr v{w+w) . (3) 

The parameter iV can be considered as the dimensionless volume of spacetime; v is the 
potential function. It is reasonable to assume the spacetime freedoms are incorporated to 
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the complex eigenvalues of the random matrix W, because the Dirac operator 

/ iW \ 

d=\ + 

^ iW + J 

invovles partial derivatives on spacetime coordinates. This is the essential assumption of the 
present work. Integrating out the gauge freedoms, we get the partition function in terms of 
eigenvalues 

r N 

Z B (N)= / dAljA^E^JJlAi-A/, (4) 

"* i=l i<j 

where A,s are complex eigenvalues(see chapter 15, [3]). It has been shown that the micro- 
scopic eigenvalue correlation at (which is the normalized correlation such that the mean 
spaceing is unit at the point) is universal. So we set v = Tr(WW + ). 

IV. ASYMPTOTIC FREEDOM PROPERTY 

The asymptotic behavior of the coupling constant in strong interaction is one of the 



most important features in QCD. In the context of renormalization group 21[, the effective 
coupling constant of QCD has asymptotic freedom, because there n/ = 6 flavors which 
makes (3(g) = —4^(11/4 — l/6n/) + 0(g 5 ). We shall provide a alternative expanation 
of asymptotic freedom using the present model. So when the energy E — > oo, when the 
gavitational effect has to arise, then expalnation is not satifactory. But in the present 
model, spacetime freedoms are the quaternion eigenvalues of the random matrix. As we 
know, as the rank of the random matrix N goes to infinitythe normalized (mean spaceing is 
normalized to be unit) two point correlation function 



fl22(r) = , ^ ( 5 ) 

7rr 



where r = I Ai — A2 1 . (see 



IQj ) • R22 goes to 1 as r goes to infinity, and goes to zero as r goes 
to zero. The bulk of this normalized two point correlation function is a constant equal to 
1. Since the fermion derterminant in the integrand of the partition function is separable for 
the eigenvalues. Therefore the eigenvalues in the bulk are statistically independent. This 
fact corresponds the asymptotic freedom in QCD. In reality, the rank N is finite, then the 
correlation functions vanish when r is large. 



V. GAUSSIAN ENSEMBLES 



A. the U(n) symmetry of Gaussian ensembles 

In the analysis of complicated nuclear spectrum, some statistical distribution on eigenval- 
ues of GOE(Gaussian Orthogonal Ensemble) are computed The eigenvalues corresponds 
to nuclear energy levels which are always positive. We shall use Gaussian complex ensembles 
to represent black holes. 

Gaussian complex ensemble is defined on the space of N x N matrices with complex 
entries by the following properties: 

(1) The probability P(H)dH that a system of Er^q will belong to the volume element 

dH =ti n dH ^ ( g ) 

7=0 l<k,j<N 

where H2,7 = 0,1,2,3 which are the four components of the quaternion Hij, is invariant 
under every automorphism 

H -> U + HU (7) 

on the space of quaternion matrices, where U is any unitary matrix. 

(2) Various linearly independent components of H are also statistically independent. 

As we know, an invariant of matrix spaces under similarity transformation only depends 
on the first N traces. However, we shall use a stonger lemma. Smooth invariants of N x 
N matrices under similarity transformation, only depends on all the N eigenvalues of a 
matrix. In fact, N x N complex quaternion matrices is a 2N 2 real dimensional manifold. 
Since nonsingular matrix space is N 2 dimensional, the action of the center of nonsingular 
matrices which are diagonal nonsingular matrices, on N x N complex matrices is invariant 
on diagonal complex matrices which is of real dimension 2N. So the quotient orbifold is 2N 
real dimensional. The orbiford is just diagonal matrix modular by the permutation group. 
We know for a continuous function defined on a matrix space which is invariant under 
similarity transformations must be a function of all the eigenvalues of the matrix. However, 
if we only assume that the similarity transformation restricted as unitary matrices, we have 
to require the function is analytic. Without proving, we state the following lemma, p 
is a analytic function on the space of complex matrices, which is invariant under unitary 
similarity transformation. Then p only depends on the eigenvalues of a complex matrix. 
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We generalize Porter and Rosenzweig theorem P(H) is an ensemble on the space 

of N x N matrices with complex entries, then P(H) is a Gaussian complex ensemble, if and 
only if 



P(H) = exp(-a tr{HH + ) + bRe tr(H) + c) 



(8) 



where a, b, c are real and a is positive, c is a normalization constant 

In fact, since P is invariant under unitary transformations, we shall consider how H is 
changed under rotation transformations which is a subgroup of unitary matrices. Let 



where 



U 



\ 

Differentiating the above equation, we get 

8H _ 8U+ 
~89 = " 



H = U + H U, 

( cos 8 sin 6 

— sin 9 cos 9 

1 





(9) 



o\ 




1/ 

,dU 



(10) 



69 HU + U+H 09 



(11) 
(12) 



And 



T dU 

89 [ 89' 



( -1 \ 
10 



(13) 



\ / 



Since each entry is independent random variables, 



p (^)=n n mm 

7 l<ij<iV 



(14) 



By the definition of Gaussian complex ensemble, P{H) is invariant under the transform of 
U, which yields 

dh\P(H) 



89 



89 



0. 



(15) 
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If we write it explicitly, we have the following equation 

+ £ D-^^ + + £ £(-^^ 2 + ( 16 ) 

7 fc>3 7 fc>3 ° kl ° k2 

These for sum terms depend on different variables. So each of them is a constant and the 
sum of the constants is zero. For each 7, we have an equation 



dH\ 2k dH. 



HI + -~7fFm k = C, (17) 

l lk ° n 2k 

Where C is a constant. C must be zero because the first term of (2.11) is a product of H] k 
with a function of Hj k , and the second term is a product of i?7i£ with a function of H^. 
Therefore we can write 

f lk = eM-a(Hl k ) 2 + cUV, (18) 

where a is a positive real constant and c(f^ k ) is a real constant. By the transform of switching 
the ith and the jth columns and rows (i < j), we can get that 

^. = exp(-a(^) 2 + c(^.)), (19) 
and since the sum of the first two terms is zero, we have 

/2 = exp(-afl2 + 6fl3 + c(/2)), (20) 

where b and c(f^) are constants. 

By fixing the unitary matrix U to be a diagonal matrix with the entry at the ith column 
and ith row equal to i, j, or k, we get for any i < j, 7, up to a constant factor 

^. = exp(-a^), (21) 

and 

^ = exp(-a(i^) 2 + ^), 
& = exp(-a(/4) 2 ), 
fix = e w{-a{H?i), 

ft = exp(-a(^). (22) 
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Therefore we are able to write P{H) by 

P(H) = exp(-a trHH + + b Re(trH) + c), 



(23) 



where c is a constant. 

Conversely, if P{H) = exp(— Na tr(HH + ) + Nb Re(trH) + c), it is straight forward to 
compute the correlation function of any two components of two entries or one entry. Then the 
computation shows that all the components of all the entries are statistically independent. 



B. joint distribution of eigenvalues 

The n point correlation function is defined by Ql( 

N\ f 

R{Xi, A 2 , • ■ • , A n ) = 7^ — J Pn(Xi, ■ ■ ■ , \^)d\ n+ i ■ ■ ■ \ N . (24) 



As shown in [10||, the normalized two point correlation fuction(in which the mean spacing at 
the origin is normalized to 1), is equal to 1 — , Here we can prove for n point correlation 
fuction goes to a constant, if |r$ — Tj\ — > oo for all it goes to zero, if |rj for all 

i,j. In fact, the normalized n point correlation function 

JV-l 

(2N) 1 / 2 



Rnorm(Xl,X 2 , ' • • , X n ) = det( T -— ^ 4>k{^i)4>k{^j))l<i,i<n (25) 



fc=0 



As iV oo, * /2 J2k=i 0fc(^i)0fc(^i) -> ^^J^, where r = (2JV)1/2 |x f — xcj*|- Therefore we 



(2iV) 

get the conclusion. 



C. how the U(N) symmetry is breaking 

If the matrix elements are exactly gaussian, as we have proved, then the gaussian complex 
ensemble has a U(N) symmetry. However, physically the matrix elements are not necessarily 
gaussian and this fact doesn't change the correlation function(Wigner's conjecture). When 
all the eigenvalues are small in a ball of origin, U(N) symmetry is approximated; when the 
eigenvalues are far from the origin, the U(N) symmetry is then broken. 
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VI. OBSERVABLES OF BLACK HOLES 



In this model, gravitational phenomena can be understood as following: we consider two 
gravitational systems in one Gaussian unitary ensemble, their spacetime distance observables 
correspond to the norm of the substraction of eigenvalues of this ensemble. When the 
eigenvalues' norm are sufficiently large, the joint distribution probability becomes smaller, 
then the classical gravitational phenomena arises; If the eigenvalues are restricted in a small 
vicinity of the origins, then because of the existence of the fermion determinant, there is a 
strong repulsion of the these eigenvalues from 0, which corresponds the black hole radiation. 

We are going to find the thermodynamics of black hole physics in terms of random matrix 
theory. Classically, there are four thermodynamical laws in classical black hole physics which 
clearly indicate the thermodynamical nature of quantum gravity. A static neutral black holes 
is determined by its mass M, angular momentum J and angular velocity Qh we have the 
generalized Smarr formula 



where A is the area of the horizon. In our present model, the symmetry is U(N) which 
is maximized. Therefore the angular velocity is zero. Since in the present random matrix 
model, all the matrix elements are statistically independent, then we know the entropy 



where Ci is the constant which equals to the entropy of each matrix element. In classical 
black hole thermodynamcis, 

< 28 > 

where r is the radius of the black hole. 

In information theory, for smooth probability distribution, there are alternative definition 
of entropy which is called relative entropy or KullbackCLeibler divergence. But here we use 
another definition because it is approximate to the shannon entropy when a is not large. 




(26) 




(27) 
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The continuous version of shannon's entropy can be negative. In fact, 
I(P(H)) = - J P{H) \nP(H)dH 

= - J exp(-atr(M+) + bRetr(H))(^-) 2N2 exp 



,Nb 2 



4a 



Nb 2 

-atr(HH + ) + bRetr(H) + 2N 2 In a - 2N 2 In vr + ) f [ J] dH? 



4 a 1 11 11 

7=0 l<i,j<N 



. _L -2iV 2 ln- 

2 4/V2 ^ 



In the following, we assume a = 1, 6 = 0. 

In the present model, if we define the mass M and r 2 proportional to 

=NC N2 I \ Zl \ 2 e~^^\ 2 Hlzi-ztfJlQzi-ztflzi-ztfdz, (29) 

^ C i i<j 

where C^2 is the normalization constant of GUE. If we put the parameter a in front of 
Y2iLi \ z i\ 2 m t ne joint probability distribution which is denoted by P(a), then its normaliza- 
tion constant is proportional to a Ar2+7V A^!(27r) iV Y[f r(2j)(see chapter 17,[l0|), i.e., 

N 

(P(a)) = const a N2+N N\{27c) N Y[T{2j). (30) 

i 

Here the constant doesn't depend on a. Differentiating both sides, we get 

M ~ r 2 ~ iV 2 (31) 



A. holographic principle 

Gerard 't Hooft proposed that in quantum gravity theoryl, the information of a gravita- 
tional system in spacetime is totally determined by the area of the surface [6 ! ]. His conjecture 
is based on the entroply analysis of classical black hole theromodynamics. From the present 
model, the information or the degrees of the freedom of the underlying system is equal to 
the product of iV 2 with a constant, which is proportional to the surface area A = 47rr 2 , 
where r here is the schwartz radius in. 
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B. the instability of Kerr black holes 



In this stationary black hole model, we have seen that the black holes are neutral, which 
means a black hole will evaporate its charge first than other matters. This is because if a 
charged black hole is in a stationary state, it has SU (N) symmetry. However the proof of 
the theorem indicates that the SU (N) symmetry enhanced to a U (N) symmetry. Thus this 
black hole in fact is neutral. 

VII. CONCLUSION 

The issue of black hole is dealt with in the context of random matrix theory in this paper. 
Spacetime geometry and degrees of freedom are combined together by the virtue of Gaussian 
complex ensembles. This model also is also an asymptotically free theory. Classical black 
hole observables can be computed in this model and holographic principle is also preserved. 
However, how the four dimensional spacetime freedom arise is related to the U(n) symmetry 
breaking process is an important issue. We would like to expect and continue to work on 
this problem. 
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